Is there more than one thermal source? 
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Abstract. BRAHMS has the ability to study relativistic heavy ion collisions over a 
C\| ■ wide range of pr and rapidity . This allows us to test whether thermal models can 

O \ be generalized to describe the rapidity dependence of particle ratios. This appears to 

work with the baryo-chemical potential changing more rapidly than the temperature. 
Using fits to BRAHMS data for the 5% most central Au+Au collisions we are able 
to describe S and Q ratios from other experiments. This paper is dedicated to Julia 
Thompson who worked to bring South African teachers into physics. 
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1. Introduction 

The purpose of RHIC is to map the phase structure of QCD. So far the community has 
concentrated on AuAu, dAu and pp collisions at ^/snn = 200 GeV in the hope of finding 
the quark gluon plasma. BRAHMS' special contribution has been to study the hadrons 
produced in these collisions over a broad range of pr and rapidity The distribution 
of particles in rapidity and px may give information on the transverse and longitudinal 
flow while the mix of different kinds of particles may tell us about the "quark chemistry" 
of the system. It is on this aspect of the bulk behavior that we shall concentrate in this 
^ ■ paper. 

2. Particle Yields and Transverse Momenta 

BRAHMS has measured particle spectra over a very wide range of rapidity and pr- Our 
AuAu spectra are summarized in Fig.Q which shows the rapidity densities, dN/dy, and 
the mean transverse momenta, (pt) , for 7r, K, p and p as a function of rapidity [2] . For 
7T, k and p the yields peak at y=0 and drop significantly at higher rapidities. The tt + 
and 7T~ yields are nearly equal while an excess of K + over K~ is observed that increases 
with rapidity. 

The mean pts, (pr) , for all particles decline slowly with rapidity. We have fitted 
the spectra themselves to the "blast wave" model in order to estimate the transverse 
flow and kinetic freeze-out temperature [3 j . We find that the flow decreases with rapidity 
while the kinetic temperature tends to increase [3]. 
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Figure 1. Rapidity densities 
(top) and mean transverse 
momentum (bottom) as a 
function of rapidity. The lines 
show Gaussian fits. 



3. Rapidity Dependence of Chemical Freeze-out 

Traditionally thermal analyses have assumed that there is only one source in heavy ion 
collisions and have used particle yields integrated over all phase space as input. At 
RHIC several groups have shown that such models can give an excellent description 
of a large number of particle ratios at mid-rapidity The use of mid-rapidity ratios has 
been justified on the grounds that there is a boost invariant region near y=0. Such 
a hypothesis implies that there is a different source away from mid-rapidity. We have 
investigated this idea by fitting our particle yields at several different rapidities to a 
chemical model [Sj. 

Figure^lshows that as y increases both the baryo-chemical potential and (to a lesser 
extent) the chemical freeze-out temperature increase. This may suggest that the system 
has fewer degrees of freedom at higher rapidities. Figure 121 suggests that chemical freeze- 
out temperature changes more slowly with rapidity than the baryo-chemical potential. 
It is possible that this is also true as we change the collision energy? Figure 01 shows 
the correlation between K~ / K + and p/p ratios measured by BRAHMS for AuAu and 
pp collisions at several rapidities. Lower energy results are also shown. Both the AuAu 
and pp results can be described by a power law, K~ / K + = {p/p) a with a = 0.32 ± 0.4 
for pp and 0.28 ± 0.06 for AuAu Ej. Preliminary AuAu data from ^/snn = 63 GeV 
are also consistent with this fit jjj. 

Recasting the power law fit in terms of chemical potentials gives /i s = (0.28±0.6) -\i q 
independent of temperature. If we take this seriously then we should be able to predict 
the correlation between other strange particle ratios and p/p = e~ 6 ^ q ^ T just by counting 
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Figure 2. Preliminary thermal 
(right) fits to BRAHMS AuAu 
data at various rapidities. The 
curves are the one sigma con- 
tours. 
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Figure 3. Correlation between 
K~ /K + and p/p ratios for cen- 
tral AuAu and pp collisions. The 
NA44 and NA49 the results are 
for central PbPb collisions. 



the number of strange and non-strange quarks. For example 
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Figure 4. Cl/Q and H/S ratios from 
STAR and NA57 versus predictions Eqns 
[1] and [2] E] The 63 and 200 GeV 
data are still preliminary. The p/p ratios 
have been displaced slightly for clarity. 
The bands represent the error on the 
determination of /i s from the BRAHMS 
data. 



and 

,=, / _\ 0.52 

-(4 M4 +2m b )/T _ „-(1.12+2)./ lg /T _ [ P \ /g) 



Figure 0] shows the predictions of Equations (1) and (2) as well as preliminary S and 
ratios from NA57 (v^^ =9 and 17GeV) and STAR (y/S^ =63, 130 and 200GeV) 
[HI E] • The agreement between the prediction and data is very good. Note that our 
ignorance of the temperature T does not affect the shape of these curves. 



4. Summary and Conclusions 

Thermal descriptions are a powerful way to describe AuAu collisions. Chemical analysis 
of our particle yields hint that both the baryo-chemical potential and the chemical freeze- 
out temperature increase with rapidity. One could interpret this in terms of the system 
becoming less partonic (with consequently fewer degrees of freedom) at higher rapidities. 
At this stage however the data are also consistent with a constant temperature at 
different rapidities. The change of baryo-chemical potential seems to dominate the 
rapidity (and energy) dependence of the particle ratios. However the correlation between 
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K~ / K + and p/p is different for pp and heavy ion collisions. This may be an effect of 
the small size of the pp system. For central AuAu (PbPb) collisions the relationship 
fi s = (0.28 ± 0.6) • fj, q derived from BRAHMS data gives a good description of hyperon 
ratios. 
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